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<N ; The pressure of a system in thermal equilibrium is expressed as a mass inte- 

gral over a sum of thermal propagators. This allows a Dyson resummation 
J> . and is used to demonstrate that potential infrared divergences are rendered 

ON ' harmless. 

oo : 

The perturbative treatment of finite-temperature quantum field theories involving massless bosonic 
fields always faces infrared divergences starting at a certain loop order. In the case of the partition 
function or pressure, these set in at three-loop order in four-dimensional theories. One relatively 
4h ■ harmless class of these infrared divergences has to do with the appearance of thermal masses. Re- 
peated self-energy insertions in a closed loop give rise to arbitrarily high powers of products of these 
thermal masses and massless propagators all with the same momentum. However, once these thermal 
masses are resummed into the propagators, such infrared divergences are avoided. As a remnant, the 
resummed perturbative series is found to be nonanalytic in the conventional loop expansion parameter 
- instead of being organized in, say, powers of coupling constant squared, the series then involves 
single powers and logarithms of the coupling. 



On the other hand, if there are modes for which no thermal masses are generated, there is a different 
source of infrared divergences at and beyond four- loop order in renormalizable theories. This is the 
case in nonabelian gauge theories at finite temperature, where perturbation theory fails to produce 
screening masses for static chromomagnetic fields. Since in nonabelian gauge theories there are vertices 
connecting exclusively magnetostatic modes, it is generally assume that perturbation theory has 

to stop at three loop order. Perturbation theory seems to be capable only of establishing the matching 
to the nonperturbative sector, which is that of an effective lower-dimensional gauge theory and which 
has to be treated by nonperturbative methods, for example on a latticed . 

An open problem of nonabelian gauge theories at finite temperature is therefore whether something 
like the generation of a thermal screening mass is indeed the mechanism which removes the infrared 
divergences of perturbation theory. But in this note we show that, in the case of the pressure, 
a screening mass is not essential to make the infrared divergences harmless. (Within the ladder 
approximation, something similar has been put forward as a conjecture in Ref. 4.) 
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The pressure at temperature T is given by 



P=^logZ (la) 



where 

logZ= 1/2 ( 1+1/12 { j +1/8 ( I ) +1/48 fcT 1 +1/8 








+1/24 — 1 — +1/16 +1/16 +1/8 



(16) 

By naive powercountingW , all the diagrams appear to be dangerous beyond three-loop order, and some 
of them already at three-loop order (if there are unresummed thermal masses). Since they all can be 
looked at as some subdiagram(s) inserted in one propagator closed onto itself, one might expect that 
a Dyson resummation could take them into account. But this expectation is spoiled by mismatching 
combinatorial factors: making a single self-energy insertion in the first diagram of log Z gives 



1/2 



(*) - w © +,/4 CO +v12 




+1/4 — H +1/8 




(lc) 

In the traditional ring-resummation^] method, one considers multiple self-energy insertions. This 
does give the correct combinatorial factors, but only for those diagrams having at least two self-energy 
insertions, which is sufficient for eliminating the infrared divergences caused by unresummed thermal 
mass insertions. The case of only one irreducible self-energy insertion is not covered and so we need 
some other method to see that there is no infrared problem when there are propagators that are not 
sufficiently screened. 

In the following, we shall derive an unconventional formula for the partition function of a quantum 
field theory at finite temperature, which has the form of a resummed one-loop quantity, while being 
exact. 

The partition function Z may be written! 2 ! as a path integral: 

Z(T) = f n#rexp(iS[0,T]) (2) 

Here, the thermal action is written in terms of unrenormalised fields and couplings: 

?] = j c (\ X>0r-d</v - m^l) + C INT ^ (3) 

In the variant of the real-time thermal field theory we are using, the integration is over all x and over 
the time-contour C which, in the complex t-plane, runs along the real axis from — oo to +oo, back to 
— oo, and then down to — oo — i/T. The action has been written as if all the fields 4> were real scalars; 
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the modification to the actual case of QCD, with quarks, gluons and ghosts, provides no essential 
difference to what we are going to do. We do need initially to include a mass for each field; in the 
end, all the masses will be set at their physical values. 

Simply giving masses to gauge fields, even only temporarily, might appear to be a dangerous proce- 
dure. Explicit mass terms destroy the gauge invariance. If we want to preserve unitarity, the ghost 
contribution changes discontinuously when masses are introduced: the massless limit is no longer 
smooth^]. However, we shall assume that the gauge theory has been gauge fixed first and the only 
change we shall make is to add masses to all the fields. The mass terms then do not change the total 
number of degrees of freedom. It does not matter that, in the presence of the masses, the theory is 
not physical; the masses are just a mathematical device. We discuss this further below. 

A simple differentiation with respect to the auxiliary masses gives 



Here, < . . . > denotes a thermal average. Space-time translation invariance tells us that the thermal 
average of 4%(x) is independent of x, and so the x integration is trivial: 




(4a) 





(46) 





(4c) 



Here 




(5) 



Dl 2 (q,T) is an element of the familiar 2x2 matrix propagator D r (g, T) of Keldysh-contourt 2 ]'^] 
real-time thermal field theory. This matrix propagator has the general structure $]M 




(6a) 



where 



D r (q,T) 



1 




q 2- m 2 r -U r (q,T) 

/ e M/2T e -w 



(66) 



with n((jj,T) the Bose distribution 



n(uj,T) = 



1 



(6c) 



e M/T - l 



On multiplying out the matrices in (6) we obtain 



Dl 2 {q,T) = -2n(q°,T) e \i°\/^ e ~ q ° /2T Im 



(7) 
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Insertion of (7) into (4) gives, with (1), 



dm 2 . 



op 



-J ■§^ N ^°^ I ^ T ) 




I r {q,T) = -Im 



1 



(8) 



q 



,2_ m 2_n r ( g ,r) 



Notice that, although we have not written it explicitly, both Tl r (q,T) and therefore I r (q,T) depend 
on all the masses. We assume that the pressure vanishes if all the masses are set to oo, because then 
all the particles are static. In order to complete the calculation of the pressure, we may set all the 
masses m r equal to the same value m, and integrate. Then the true pressure, for all the masses zero, 
is 



We have checked explicitly, in the first few orders, that (9) does reproduce the diagrams in (lb) with the 
correct combinatorial factors. This can be seen by expanding out the Dyson resummed propagators, 
upon which the mass integral can be eliminated by appropriate partial integrations. This also works 
for the more precariously looking case of nonabelian gauge theories, where Feynman gauge turns out 
to be the simplest case to consider. That we are spoiling BRS invariance by giving masses to all 
fields does not appear to be a problem. This should not come as a surprise, since the diagrammatic 
expansion in (lb) does not know about BRS. However, we clearly need to have a gauge fixed theory; 
otherwise even the interaction- free pressure would come out wrong. M 

We now argue that (9) is not infrared divergent. First, it is evident that any infrared divergences 
that may be present in the self energies II r ((/, T, m) for some values of its arguments will not cause a 
divergence in P; the Dyson resummation, which puts U r into the denominator of I r , has seen to this. 
Rather, what might cause trouble is a zero of the denominator of an I r , or of N(q°,T). 

What is responsible for the increased sensitivity to the infrared regime at finite temperature is the 
singularity of N(q°,T) at q° = 0. Near q° = 0, it behaves as e(q°)T/q°. However, this is rendered 
harmless by I r . e(q°)I r (q,T) is proportional to the spectral density defined by 






(q,T). In fact, the integral over 



(10) 




(11) 



Using this, we can write (9) as 



f'OO f' 



d 3 q 1 



(2vr) 3 q 2 + m 2 + U r (q° = 0, q, T) 
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+E/ °^ 2 /(£y T ( Ar(9 °' T) "^) Ir[q ^ m) (12) 

which separates off what would be identified as the zero-mode contribution in the imaginary-time 
formalism. 

The zero-mode contribution is clearly harmless in the infrared, i.e. |q| — > 0, as long as the number 
of spatial dimensions is greater than two. This is true even for m = and for the case where II r 
does not furnish a mass term. At q° = 0, one would not expect singularities at real values of |q|. 
Nevertheless, such singularities can appear at some stages of a perturbative evaluation of II r . In 
the magnetostatic sector of QCD, these spacelike singularities have been termed "Landau ghosts" 

[10] 

These are generally assumed to disappear in more accurate calculations of IT r , but even so their 
presence would not introduce divergences in (12). The corresponding singularities are smeared out by 
the additional integration over m 2 . 

Actually, the zero- mode contribution itself is ultraviolet divergent, but this divergence can be avoided 
by using dimensional regularization^] ; otherwise it has to be combined with the second term. Fur- 
thermore, there is a logarithmic divergence in the integral over m 2 for large m 2 , but this is merely an 
artifact of the above separation. In the complete expression (9), I r behaves as 1/m 4 when m 2 — > oo, 
assuming that Imll,. does not blow up in this limit. 

Consider now the singularities that might arise in (12) from a zero in the denominator of an I r . Clearly, 
the denominator can vanish only if H r (q,T,m) is real. This will be the case for q° = 0, whereas for 
nonvanishing qo it will happen only in some low order of a perturbative evaluation. Even then, II r 
produces a mass gap (the plasma frequency), which protects us from infrared divergences also when 
m = 0, while m > will only improve the infrared behaviour. Moreover, when evaluated through 
higher orders, n r will eventually produce an imaginary part, which removes the singularity from the 
real axis. 

Since, in the form (9), the pressure is manifestly infrared finite, this might provide a better starting 
point for a correspondingly improved perturbative evaluation of the pressure than standard pertur- 
bation theory, which, in nonabelian gauge theories, is expected to have an insurmountable barrier 
beyond three loop order. In (9), the primary objects to compute are self energy diagrams. As long 
as these are kept in the denominator of (8), all infrared problems should be eliminated. When this 
becomes essential, one clearly is going beyond ordinary perturbation theory: the final result will not 
have a power-series expansion in the coupling. Of course, a practical evaluation of (12) will have to 
cope with a great number of technical difficulties, but it does seem to offer a way to circumvent the 
no-go argument against a perturbative approach to the pressure. 
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